Generating dark solitons by single photons 
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We show that dark solitons in ID bose systems may be excited by resonant absorption of single 
quanta of an external ac field. The field frequency to should be slightly blue-detuned from e s (hq)/h~, 
where e 3 (hq) is the energy of a soliton with momentum corresponding to the external field wavenum- 
ber q. We calculate the absorption cross-section and show that it has power-law dependence on the 
frequency detuning u) — e s (hq)/h. This reflects the quantum nature of the absorption process and 
the orthogonality catastrophe phenomenon associated with it. 
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The existence of dark solitons (DS) is probably the 
most spectacular manifestation of the role played by weak 
inter-particle interactions in ID cold atomic gases [lj]. 
Such solitons are macros comically large areas of partially, 
or even completely depleted gas, which propagate with- 
out dispersion. It is natural to interpret these objects as 
solutions of the classical Gross-Pitaevskii equation 0, 0| . 
Correspondingly, the means to create DS, employed so 
far, required a macroscopic classical perturbation applied 
to the atomic cloud. An example of the latter is the phase 
imprinting technique where a finite fraction of the 
ID cloud is subject to an external potential for a certain 
time. 

Yet, there is a deep connection between DS and intrin- 
sically quantum nature of the gas. It was first appreciated 
by Kulish, Manakov and Faddeev [f|, who noticed that 
the soliton's dispersion relation coincides @ with the ex- 
act lower bound of the quantum many-body spectrum of 
the Lieb-Liniger model [7], Fig.Q] The existence of such a 
bound, known also as the Lieb II mode, was derived ear- 
lier using Bethe Ansatz technique Q. This observation 
means that the DS may be considered as a classical ap- 
proximation of very peculiar quantum many-body eigen- 
states: those which possess the minimal possible energy 
62 at a given momentum hq, see Fig. [TJ 

In this Letter we show that the underlying quantum 
nature of the gas allows for a qualitatively different way 
of creating DS. Namely, the soliton may be generated 
by the absorption of a single quantum |8| of an ac field 
slightly blue-detuned from the resonance with the soliton 
energy fko > £2(^9), where lo and q are frequency and 
wavenumber of the ac field. Notice, that the comparison 
between the soliton's energy and the ac frequency does 
not appear in the classical treatment at all. 

The quantum efficiency of the discussed process is lim- 
ited by the orthogonality catastrophe [9| . Indeed, a state 
of the system after it is excited from the ground state 
by absorbtion of one quantum of the ac field has a very 
small overlap with the eigenstates forming the soliton. As 
a result, the quantum yield exhibits the very characteris- 
tic power-law dependence on the energy excess above the 
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FIG. 1: Energy vs momentum plane for excitations. A photon 
is represented by a solid dot. The "light cone" with Bogoli- 
ubov velocity vb determines the range of possible momenta 
of generated dark solitons, p_ < p s < p+. 



threshold: ~ (tko — £2)^ ■ The corresponding exponent 
fi2 is a function of the wavenumber q and the dimen- 
sionless interaction parameter 7. For weak interactions, 
7 <C 1, the exponent is large, [i 2 ~ 7 -1 / 2 ^> 1, reflecting 
the strong orthogonality. The photon absorption and DS 
production, associated with it, are higher for a moder- 
ate interaction parameter, 7^1, when the DS density 
depletion is weaker. 

To set the notation let us briefly remind what the local- 
ized solutions of the non- linear Gross-Pitaevskii equation 
are. Quasiclassically, the condensate wave function obeys 
the following equation [Io| : 



id t ^ + — a 2 * + c in - |*| 2 ) * = , 
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where n = N/ L is the average concentration and L is the 
length of the system with the periodic boundary con- 
ditions. The interaction strength c gives rise to the 
dimensionless parameter 7 = cm/ (fin) whose smallness 
7 <C 1 is the criterion of the weakly interacting gas. 
Soliton solutions of Eq. |T]) have the form 11 1 
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where the only free parameter 9 S is the change of phase of 
the condensate wave function = y/n(x)e t6 ^ across 
the soliton. The condensate density n(x) reaches equi- 
librium value n away from the soliton (see Fig. [2]). For 
a given 9 S the soliton velocity and its spatial extent 
arc fixed and given by v s — vb cos(9 s /2) and l s — 
H (itivb sin(# s /2)) _ , where vb denotes the Bogoliubov 
sound velocity, which in the weakly interacting gas is 
given by vb = y^nh/m. The soliton, Eq. |2]), repre- 
sents the density depletion, propagating without disper- 
sion with the velocity v s . The number of particles pushed 
away from the soliton core is 
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where K = irnh/ (mvs) is the thermodynamic compress- 
ibility of the gas. 
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FIG. 2: (color online) Density n(x) and phase -&(x) profiles of 
a soliton in a system of length L, as given by Eq. @. Note 
that the density perturbation is local, while the perturbation 
of phase is not. 

In the presence of a soliton the system acquires a non- 
zero momentum p s and energy e s which are expressed 
through the phase shift 9 S as: 



p s = hni 



sin 9 S ) 



-finvssin 3 y . (4) 



These two relations implicitly define the soliton disper- 
sion relation e s = e s (p s ). One may check that the ve- 
locity of the soliton satisfies v s = de s /dp s , as expected 
for a particle. It is important to stress that while the 
energy e s is entirely associated with the soliton core, the 
momentum p s is shared between the core and the rest of 
the gas. More precisely, the soliton core carries the mo- 
mentum — hn sin 9 s (the negative sign here corresponds 
to a hole in the density), whereas the momentum Hn9 s 
is uniformly spread to the rest of the gas. The latter 
fact is easily seen from the phase ^(x) profile, Fig. 
and the observation that the momentum density of the 
condensate is given by hnd x $. 

At small momenta p s -C nhn, one finds 9 S — * 
and thus e s » vbPs, cf. Eq. dj}. In this limit 
the soliton dispersion approaches the phonon branch 
[l2| characterized by the Bogoliubov spectrum tusiq) = 
( v Bq)\/l + {hq/2mvB) 2 , see Fig. [TJ On the other hand, 
for ps — irhn, the phase shift is 9 S = tt so the soliton is 



at rest: v s = 0. In the latter case the entire momentum 
p s is uniformly spread over the bulk of the gas. 

Imagine now that the gas is subject to a weak space 
and time dependent external potential Vo cos(gx — cut). 
According to the Fermi Golden rule, the system may ab- 
sorb quanta of this field if its many-body spectrum pos- 
sesses excited states with the momentum hq and energy 
hio. As may be learned from the exactly solvable model 
, such states form a continuum whose energy is bound 
from below by the line £2(hq) (so called Lieb II mode). 
It was subsequently noticed in Ref. [f| that in the limit 
of weak interactions 7 <C 1 the Lieb II mode coincides 
with the soliton dispersion relation, i.e. E2(hq) = e g (Hq). 
These observations imply that an absorption of a quan- 
tum blue-detuned from the soliton energy hu> > e s {hq) 
is (i) possible and (ii) leads to creation of the soliton. 
Our aim is to study the scattering crossection of such 
processes. 

To this end we notice that a photon absorption first 
creates a virtual state of the condensate with a local per- 
turbation of the condensate wavefunction. Since the pho- 
ton carries momentum hq and carries no extra particles, 
so does the initial local perturbation. Subsequently this 
perturbation evolves according to the equations of mo- 
tion and eventually it must take a form of a superposi- 
tion of real excitations, i.e. conserving overall energy huj 
in addition to the momentum Hq. We expect that such 
a final state contains a soliton with momentum p s ~ hq 
and core energy £ s {p s ) < The small excess energy 
huj^e s > and the momentum difference hq—p s are car- 
ried away by the small-amplitude sound waves (phonons) 
with velocity vb- This observation immediately implies 
that vs\hq— p s \ < (ti^—EsiPs)), with the equality reached 
if all the sound waves are emitted in one direction only. 
As a result, the range p~ < p s < P+ of possible soliton 
momenta p s is limited 



p±(q,u) ^hq± 



huj — £ s {hq) 

v B ± v„(hq) 



P+-P- < Ps 



(5) 



The corresponding construction is shown in Fig. [TJ 

We found that the inelastic scattering crossection for 
absorbing a photon with the wavenumber q and fre- 
quency u>, while creating the dark soliton with the mo- 
mentum p s is given by 
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It is characterized by the power-law dependencies on the 
deviations of the soliton momentum from the upper and 
lower kinematic boundaries p±(q, w). The corresponding 
exponents p,± are functions of the soliton parameter 9 S = 
9 s (q) and the thermodynamic compressibility K, 
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FIG. 3: (a) Inelastic scattering crossection as a function of 
the soliton momentum p a . (b) Momentum dependence of the 
exponents fi± for K = 10. The exponent in Eq. ([8)l is given 
by M2(?) = + /-*-(?) - 1- 



In the last equality here we employed Eq. (J3J), and 
the function 8 s (q) is implicitly defined by Eq. ([J]) with 
p s = Hq. The calculated crossection as a function of the 
soliton momentum p s for a photon with q = irn is plot- 
ted in Fig. [3^. We have also plotted the (^-dependence of 
exponents /i±(<z) in Fig.[3]3. 

Equations (O - ([7]) are the main result of this paper. 
They provide the probability of the soliton excitation 
with a specific momentum p s (and therefore specific ve- 
locity v s = de s /dp s ). One may also be interested in the 
integral probability of absorbing a photon with a given q 
and w > e s (hq)/h, which results in creation of a soliton 
with an unspecified momentum. This quantity is nothing 
but the dynamic structure factor (DSF) S(q, u>) of the ID 
Bose gas. Integrating W q ^{p s ), Eq. ([6]), over the soliton 
momenta p s , one finds for DSF in an immediate vicinity 
of the lower spectral boundary hw > £2(^9) 
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Being multiplied by the intensity of the radiation Vq this 
quantity gives a number of solitons excited per unit time 
and per unit length of the irradiated ID gas. According 
to Eq. (O, the wavenumber-dependent exponent ^2(9) = 
+ - 1 is given by 



2K 



(9) 



The power law behavior of DSF near the lower spectral 
boundary £ 2 (/tq) was suggested earlier in Ref. [13| . How- 
ever, Ref. [l3( derived functional dependence of ^2(9) 
only in the limit of strong interactions, 7 > 1. Equa- 
tion @ provides the answer in the opposite limit of weak 
scattering, 7 <C 1. 



Below we outline the details of derivation. As was ex- 
plained earlier, the initial local perturbation of the con- 
densate carries the momentum of the absorbed photon hq 
and no excess particles. After a short time, which may be 
estimated asr s = I s /vb, the soliton core is formed along 
with a bunch of outgoing sound waves. The core, which is 
the density depletion, carries momentum — fori sin # s and 
— N s particles. At times t > t s it propagates without 
dispersion and behaves as a free particle with the energy 
£ s . The remaining momentum h(q + n sin 8 s ) — hn9 s , 
cf. Eq. (|U), and N s particles, localized on a scale ~ l s 
at t r~j t s must be carried away and spread over the en- 
tire system at t ^> t s by the linear sound waves. The 
latter are conventionally described [3] by the linearized 
hydrodynamic Hamiltonian 



H sw = 4r^- / < 

Z7T 



\d x <t>) 2 + K{d x df 



(10) 



where hnd x d{x) and d x (pji: are operators of momen- 
tum and particle excess densities, correspondingly. Their 
canonical commutation relation reads [ip(x),'d(y)] — 
i(?r/2)sgn(a; - y). 

The injection of momentum hn6 s and N s particles into 
the sound waves at a point x is achieved by acting with 
operator 



(11) 



on the ground state of the Hamiltonian (jT0[) . Indeed, the 
operator e 1 "^^ shifts -&(y) by 9 S at x > y, accommodat- 
ing momentum Hn9 s . Similarly, e lN ^( x ) shifts <p(y) by 
N s , at x > y accommodating A^ s particles at the point x. 
Since we are dealing with large shifts, we can disregard 
non-commutativity of the operators, effectively adopting 
the semiclassical approximation (see below). 

The energy of a state created immediately upon the 
shift Eq. (fTTj) is much higher than Hcu(q) — e s (hq). In full 
analogy with the instanton picture of tunneling across 
a weak link in a Luttinger liquid [l5| . the created vir- 
tual state evolves, eventually reducing its energy to 
hu>(q) — e s (hq). To find the photon absorption crossection 
resulting in a creation of a dark soliton, one needs to eval- 
uate the time evolution of the density-density correlator 
of the excited sound waves, 



G(x,t) = ($ sw (x,t)i>l w (0,0)). 



(12) 



Its calculation uses Eqs. (fTT|) and pop and follows the 
standard route of the bosonization theory 14 1, leading 
to 



G(x,t) 



x — vst 
il s 



x + vst 
il s 
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Here the exponents are given by Eq. ([7]) and l s enters 
through the short-time cutoff t s as / s = vbt s . Finally 
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the absorption crossection is proportional to the Fourier 
transform of the correlator: 

W q , u (p a ) ~Im JdxdtG(x,t) e i^-^/H)t-( q - Pa /h)x _ 

(14) 

where we took into account the momentum p s and en- 
ergy e s (p s ) carried away by the soliton. To evaluate the 
integral in Eq. (Tj"4")) . we take into account that for small 
energy excess fko — e s (hq) <C fku the range of the al- 
lowed soliton momenta p s is rather narrow, see Fig. [1] 
and Eq. (0, and centered around the photon momen- 
tum Hq. One may therefore expand the soliton energy as 
£s(Ps) ~ e s {hq) + (p s — hq)v s . Equations ((5]) and ([6} then 
follow from the straightforward calculation. 

The developed theory is applicable in the vicinity of 
the spectrum of the Lieb II mode. The width of the 
corresponding region in the (q, to) plane, see Fig. [T] is 
determined by the condition fku — £2(^7) ^ ^/t s - Also 
u> is restricted to be below the Bogoliubov mode u>b and 
its replica [l8|, shown by the dashed line in Fig. [TJ The 
soliton production rate reaches its maximum at fku — 
S2(hq) ~ H/t s and decreases at higher frequency. 

Additional restrictions are set by the applicability of 
semiclassical approximation. The semiclassical descrip- 
tion of the Lieb II mode by dark solitons fails at p suf- 
ficiently close to or l-xTm. This already can be seen 
from a comparison [16J of the true spectrum of the Lieb 
II mode with the prediction of the semiclassical Gross- 
Pitaevskii approach. At 7 <C 1 the two spectra signifi- 
cantly differ from one another at p < Hwy 3 ^ 4 . Equations 
([7]) and @ predict monotonic decrease of the exponents 
jj,± and with the decrease of q. For q <C 7m they yield 

K f 6q \ 2 ^ 3 2 

M?) ~ M+W = — — , v-ocq <^+- (15) 
7r z \ n J 

This behavior is actually valid in the interval wy 3 / 4 <C 
q <C 7m, the lower limit here being set by the applica- 
bility of the semiclassical approximation. In this range 
the exponents /Z2,/x + 1, while there are no further 
restrictions on /i_ > 0. 

Note that in the case of strong interaction [l3[ (7 3> 1), 
the asymptote at q —* is ^2(9) oc q. In this limit DS 
degenerates into a single hole within weakly-interacting 
Fermi sea in the effective fermionic description (l7| . We 
expect the linearity must stay for the smallest q at any 
interaction strength, and thus expect a crossover from 
/i+,/i2 oc q 2 l z to the linear behavior at q < n~/ 3 / 4 . The 
similar crossover occurs in the narrow vicinity of the 
q = 27m point, where the semiclassical soliton descrip- 
tion also runs out of the applicability. 

To conclude, we have shown that DS may be gener- 
ated by absorption of a single photon of an external ac 
field. Quantum efficiency of such a process is maximized 
at slight, ~ h/T s , blue detuning of the photon energy fku 



from the soliton energy e s (hq). Within this range the 
absorption probability behaves as the power law of the 
detuning, reflecting the quantum orthogonality catastro- 
phe phenomenon. 
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